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13  ABSTRACT 


In  this  report  a  family  of  baseband  suboptimal  communication  systems^, 
referred  to  as  "the  family  of  predistorted  replica  correlation  receivers," 
is  presented  and  analyzed. 

The  resulting  receiver  structures  estimate  the  reference  waveform 
directly  from  the  received  process;  whereas,  the  well-known  matched  filter 
or  correlation  receiver  uses  a  stored  reference  waveform.  In  this  report 
the  Doppler  performance  of  two  such  structures  is  investigated,  the  results 
indicating  that  such  structures  can  be  relatively  insensitive  to  Doppler. 
The  characteristic  function  for  the  "decision  statistic"  is  computed  for 
the  discrete-time  systems  and  numerically  inverted  to  obtain  the 
cumulative  distribution  function.  Chernoff  bounds  are  computed  to 
evaluate  the  behavior  of  the  'tail  probabilities. 
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Abstract 

In  this  report  a  family  of  baseband  suboptimal  communication 
systems,  referred  to  as  "the  family  of  predistorted  replica  correlation 

receivers,"  is  presented  and  analyzed. 

The  resulting  receiver  structures  estimate  the  reference  waveform 
directly  from  the  received  process;  whereas,  the  well-known  matched 
filter  or  correlation  receiver  uses  a  stored  reference  waveform.  In 
this  report  the  Doppler  performance  of  two  such  structures  is  investi¬ 
gated,  the  results  indicating  that  such  structures  can  be  relatively 
insensitive  to  Doppler.  The  characteristic  function  for  the  "decision 
statistic"  is  computed  for  the  discrete-time  systems  and  numerically 
inverted  to  obtain  the  cumulative  distribution  function.  Chernoff 
bounds  are  computed  to  evaluate  the  behavior  of  the  "tail  probabilities. 
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I  Introduction 

In  this  report  we  investigate  the  performance  of  several  different 
communications  systems.  The  receiver  structures  for  all  of  these  systems 
belong  to  a  family  which  we  will  call  the  family  of  predistorted-replica 
correlation  receivers.  This  terminology  arises  because  of  the  differences 
between  a  "predistorted-replica  correlation  receiver"  and  the  well  known 
correlation  receiver  (or  matched  filter).  The  latter  uses  a  stored 
reference  signal  and  correlates  this  with  the  received  process,  whereas 
the  former  derives  its  "reference  waveform"  directly  from  the  received 
process. 

We  will  begin  our  discussion  by  considering  the  conventional  binary 
antipodal  pulse-amplitude  modulation  (PAM)  system,  for  which  the  trans¬ 
mitted  signal  is  assumed  to  be  of  the  form: 

00 

x(t)  »-Jm.  s(t  -  kT),  (1.1) 

k-0  K 

where  m^  e  (-1,  1}  and  s(a)  =  0  for  all  a  i  [0,T).  For  coherent 

00 

communication,  the  receiver's  task  is  to  estimate  the  sequence  (m^J^o* 

As  is  well  known,  for  the  case  that  the  received  process  is  given  oy 

r(t)  *  3  x(t)  +  n(t)  (1*2) 

where  3  is  the  channel  attenuation  constant,  and  n(t)  is  a  white 
Gaussian  noise  process,  the  optimal  receiver  is  the  correlation  receiver 
shown  in  Figure  1.1,  with  i>(t)  =  x(t),  s^^  ( (t)^)  =  s  ( (t)^)  ,  where 
(t)^,  *  t  modulo  T.  The  term  "optimal"  is  used  here  as  "the  best  that 
can  be  done  for  any  reasonable  criterion  of  goodness  without  making  use 
of  the  source  dynamics."  The  reader  is  referred  to  ([1]  -  [A])  for 
a  more  lucid  discussion  of  these  points. 

1 


— ... 
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Now,  if  the  channel  is  such  that  the  signal  component,  S(t) ,  of 
the  received  process,  r(t),  can  be  written  as 

oo 

S(t)  •  l  sd(t  -  kT),  (1.3) 

k=0 

where  s  (a)  =  0  for  all  a  i  [0,T),  it  is  easily  seen  that  the  receiver 
d 

of  Figure  1.1  retains  its  optimality  with  sref  ((t)T>,  for  the  case 
that  scj(.)  is  unknown,  one  can  either  use  sref^t^T^  *  8((t)^)  anc* 
suffer  the  resulting  degradation  in  performance,  or  resort  to  a  different 
receiver  structure,  thus  leading  us  to  the  so-called  family  of  predistorted 
replica  correlation  receivers. 
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II  Fixed  Lag  Autocorrelation  Receiver 
A.  Continuous  Time 

We  begin  our  discussion  of  predistorted  replica  correlation  receivers 
with  the  Fixed  Lag  Autocorrelation  (FLAC)  receiver.  Consider  the 
receiver  structure  illustrated  in  Figure  2.1,  which  is  a  continuous- time 
version  of  the  system  proposed  by  Farden  [5  ].  The  output  of  the  cor¬ 
relator,  y(t),  can  be  written  as 
t 

y(t)  =  /  r (t)t(t  -  T)dx 
t-T 

t  t 

-  /  5(t)S(t  -  T)dT  +  /  S(x)n(x  -  T)dx 

t-T  t-T 

t  t 

+  /  S(t  -  T)n(x)dx  +  /  n(x)n(x  -  T)dx.  (2.1) 

t-T  t-T 

If  the  signal  component,  S(t),  of  the  received  process,  r(t),  can  be 
written  as 

S(t)  =  m  s  ( (t)  ), 

[fi  d  T 

where  [x]  denotes  the  "largest  integer  contained  in  x,"  (t)T“t  modulo  T, 
sd(a)  =  0  for  all  a  ^  [0,T),  and  n(t)  is  Gaussian  White  Noise  (GWN) , 
we  can  easily  obtain 

.  t 

z(t)  «  E{y(t) }  =  /  m  m  s  ((x)  )s  ((x-T)  )dx.  (2.3) 

t-T  [^1  [^]  "  T  “  T 

The  functions  (£t)_  and  f^-1 are  illustrated  in  Figure  2.2  for 
c  e  (.75,  1.  1.25}.  Finally,  at  t  -  (k+l)T,  we  have 
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z  =  E{y(t)  }  I  =rtL  tn  (k|1)T  S  2((t)  )dT 

k  t=(k+l)T  1  kJT 


where 


=  "k  Vl  csd  <0)- 


Cs  (y)  =  /  sd(T)  sd(x-Y)dT . 


(2.4) 


(2.5) 


Carefully  note  that  C  (.)  is  not  a  cyclic  correlation  and  that 

Sd 

s  (a)  =  0  for  all  a  t  [0,T).  This  result  suggests  that  if  the  sequence 
d 

{n^}*  0  iS  3  differentially  encoded  version  of  some  information  bearing 

sequence  {b  }°°_  ,  the  FLAC  system  has  some  chance  of  success,  at  least 
k  k“l 

for  high  signal-to-noise  ratios.  The  differential  encoding  scheme  to 
be  used  at  the  transmitter  is  the  following  (for  m  *  1,  and  k =  1,2,3, ...) . 


or  conversely. 


Vi-  if  bk  *  1 


Vi-  if  bk  ■  -l- 


1,  if  ^Vl-! 


(2.6a) 


(2.6b) 


L'1* if  ™k  Vi " _1 

Now,  suppose  that  the  signal  component  of  the  received  signal  has 


suffered  a  Doppler  shift,  i.e. 


S(t)  =  m  _  s  ((Ct)  ). 

if]  d  T 

The  mean  of  the  correlator  output  (2.3)  becomes 


(2.7) 


V‘>  ■  !  "ril,  “rillza,  sd(<,*>T>8d«E-mT>d'-  (2-8) 

t-T  L  T  I  l  j  I 

Using  the  change  of  variable  u  -  ,  and  defining  x  -  ^jr  and 

z.'(x)  -  z  (t),  we  have 
^  * 
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z'  (x)  =  j  J  ,  m  s  ((u)  T)  s  ((u-O  T)du. 

C  5  x-C  11  [u-C] 

Partitioning  the  interval  [x-£,x]  as. 

a^  ■  x-£  <  a^  <  a2  <  *‘*<  aL  *  x* 

I  I 

so  that  [u]=  k+n£,  [u-C]  -  k+P£  for  all  u  e[a£  a£)  for  £  =  1,2 . L, 


where  k  =  [x]  we  may  write: 

L 


■  7  £  \N,  Vwt  J  8d«*>iT>  8<* 


Sj ((u),T)  8>,((u-£;)1T)du. 


£-1 


Furthermore,  defining  a£  ■  a£_^  -  k  -  n£,  b£  =  a£  -  k  -  n£,  and 


q£  -  n£  -  P£,  we  obtain: 


z’  (x)  =  j  J  8d(uT)8d((u+qrOT)du.  (2.9) 


It  is  easily  shown  that  for  £,  e  (1/2,  3/2),  we  have  n£  e  (-2,  -1,  0}, 

and  P  c  (-3,  -2,  -1,  0}.  Table  2.1  shows  all  of  the  parameters 
i* 

necessary  for  the  computation  of  (2.9)  for  any  known  signal  ( • )  an<* 

for  any  value  of  £  e  (1/2,  3/2). 

Whereas  in  the  absence  of  Doppler  shift,  information  was  being 
sent  and  received  at  the  rate  of  one  bit/T  seconds,  with  a  Doppler  shift 
information  is  being  received  at  the  rate  of  one  bit/^  T  seconds;  hence, 
the  receiver  must  observe  the  correlator  output  and  make  a  decision  once 
every  T/C  seconds.  Consequently,  we  will  now  investigate  when  the 
maxima  and  minima  of  (2.9)  should  occur.  Hopefully,  the  maxima  and 
minima  of  (2.9)  will  occur  for  the  values  of  x  for  which  the  quantity 


6  -  max(b?-a.) 

X  o  X,  X, 


(2.10) 


is  maximum  with  the  constraint  that  q£  Applying  this  criteria  with 


Table  2.1 


Parameter 

values  for  computing  output  of 

FLAC 

receiver 

(x^  £{•) 

£ 

(a)  1/2  <  £  <  1 

a£  b£ 

I_£ 

II 

It 

[0,  25  -  1) 

1 

(x)L  +  1-5  5 

-i 

-2 

1 

2 

5  1 

-1 

-1 

0 

3 

0  (x)L 

0 

-1 

1 

[25  -  1,  £) 

1 

(x)L  +1-5  1 

-1 

-1 

0 

2 

0  (x)L 

0 

-1 

1 

[£,  1) 

1 

(x)x  -  5  £ 

0 

-1 

1 

2 

£  (x)L 

0 

0 

m 

(x)L  e{.} 

_£ 

(b)  5  =  1 

a  b 

a£  £ 

It 

l£ 

It 

[0.  1) 

1 

(x)x  1 

-1 

-2 

i 

2 

0  (x^ 

0 

-1 

i 

(x^  e{.) 

£ 

fc)  1  <  5  <  3/2 

a  b 

a£  £ 

It 

l£ 

It 

o 

«• 

1 

h-» 

1 

(x)  +2-5  1 

-2 

-3 

i 

2 

0  5-1 

-1 

-3 

2 

3 

£  -  1  1 

-1 

-2 

1 

A 

o  M1 

0 

-2 

2 

Ic-l.  2C-2) 

1 

(x).  +1-5  £  -  1 

-1 

-3 

2 

2 

£  -  1  1 

-1 

-2 

1 

3 

0  5-1 

0 

_2 

2 

A 

£  -  1  (*>! 

0 

-1 

1 

[25  -  2,  1) 

1 

(x)  +1-5  1 

-1 

-2 

1 

2 

0  5-1 

0 

-2 

2 

1 

2 
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the  aid  of  Table  2.1,  we  find  that  (i)  for  1/2  <  £  <  1,  (3^  attains 
its  maximum  value  of  £  for  (x)^  *  £,  (ii)  for  £  "1,  8^  attains  its 

maximum  value  of  1  for  (x),  =  0,  and  (iii)  for  1  <  £  <  3/2,  8 

attains  its  maximum  value  of  2  -  £  whenever  (x)^  e  [0,  2£  -  2]. 

Using  these  results,  we  find  that  for  1/2  <  £  <  1  and  x  *  k  +  £: 


z'  (x)  -  j  _  1  I  sd(uT)sd((u  +  i  -  £)T)du 


£  “k  mk  -  1  Csd 


(2.11) 


where  the  last  equality  follows  from  (2.5)  and  the  fact  that  sd (a)  =  0 
for  all  a  i  [0,  T) .  Similarly,  we  find  that  for  1  <  £  <  3/2  and 
(x)x  e  [0,  £  -  1): 

T  \ 

z'  (x)  =  7  m.  9  m  _  /  s  (uT)s  ((u  +  1  -  £)T)du 

£  5  k  -  J(x)i  +  2-^  d 

Cs  <(5-2,T)  +I"k-l,\-2  C,  <<5-1)T) 


+  *£  “k  2  I  sd(uT)sd((u  +  2  -  £)T)du,  (2.12) 

whereas,  for  1  <  £  <  3/2  and  (x)^  e  [£  -  1,  2£  -•  2],  we  obtain: 


sd(uT)sd((u  +  2  -  £)T)du 


ZC  (x)  "  7  ^-l  “k  -  3  /  s  (uT)s  ((u  +  2-£)T)du 

t.  £  K  1  J(x)1  +  l-£ 

+  H-1  “k-2  C.«?*1)T)+i\V2  C  ((£-2)T) 


+  ■-  m^  m^  _  ^  /  sd(uT)sd((u  +  1  -  £)T)du.  (2.13) 

Finally,  a  worst  case  analysis  of  both  (2.12)  and  (2.13)  (m^  _  ^  *\_2 
\-2Bk-3‘\-l\-3"\\-21'1  in  (2>12)  and  \-l  \-2  "  ^ 


\-lmk-3"mk%-2  =  mk  Vl 


-1  in  (2.13))  yields 


zi  (x)  >  7  Co  ( (£  -  1)T)  -7  { (£  -  1)P  +  |C  ((£  -  2)T)|},  (2.14) 

r  —  r  a  r  max  S  «■ 
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where  we  have  defined 


P  „  “  T  •  max  s  (uT) , 
max  u  e [0, 1)  d 


so  that  for  b  -  a  <  1,  we  have 


/  ■1(»T)..«u  -  Y)T)du|  P  . 


max 


G  *  7  |C  ((€  -  1)T)|/C  (0), 

^  Sd  Sd 


Now  define  the  gain  of  the  system  against  Doppler  to  be  G  -  |z'  (x)|/C  (0) . 

C  8d 

From  (2.11)  with  1/2  <  £  <  1,  x  =  k  +  £  and  *  1,  we  have 

(2.15) 

“d 

From  (2.14),  for  1  <  5  <  3/2: 

0  i  FcSoT  <  lcs  <«  -  »T> I  *  « -  »>'««  -  lcs ,  <«  -  i > 

s,  d  d 

(2.16) 

Suppose  for  a  first  approximation  that  s_(.)  has  time  bandwidth 

d 

product  TW  and  that  C  (y)  is  given  by: 


C  (Y) 

Sd 


\CB  (°)  (1  -  K|W).  for  |  Y  |  <  ^ 


,  elsewhere 


(2.17) 


For  this  choice  of  C  (.),  (2.16)  becomes,  for 
l  Sd 

and  1  <  C  £  1  +  — 

C  >  3-~  ^  -  1)  TW  . 

K 

whereas  (2.15)  becomes,  for  1  -  <  £  <  1, 


P 

max 


2C  (0),  TW  >  2, 
8d 


(2.18) 


r  „  1  -  (1  -  g)TW 

€ 

Figure  2.3  illustrates  equations  (2.18)  and  (2.19)  for 
20,  50}. 


(2.19) 
TW  e  (2,  5,  10, 
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B.  The  Discrete  Time  FLAC  System 

We  will  now  discuss  a  discrete  time  version  of  the  FLAC  system. 

The  carrier  signal  s(t)  will  be  extended  for  all  t  >_  0  via  the  rela¬ 
tionship  s(t)  *  s((t)T).  For  definiteness  s(t)  will  be  considered 
to  be  a  sampled  and  held  version  of  a  pseudorandom  noise  sequence  (PRS) 
which  assumes  the  binary  values  of  +1.  Such  sequences  are  easily 
generated  by  simple  shift  register  circuts,  and  have  desirable  auto¬ 
correlation  characteristics  [6].  The  value  of  s(t)  for 
kNT  T  <  (k  +  1)  NT,  k  *  0,  1,  2,...,  will  be  denoted  as  s^,  where 
N  and  T  are  positive  numbers  whose  significance  will  become  apparent 
momentarily.  With  this  notation,  the  autocorrelation  of  s  is  defined 


as 


0,  i  j1  nK,  n  *  0,  +1,  +2,... 
K,  i  =  nK 


(2.20) 


where  K  is  the  period  of  the  PRS,  i.e.,  s^  -  s^+K  ^or  ^  —  0. 

The  actual  value  of  K  is  determined  by  the  bit  length,  BL,  of  the 

BL 

shift  register  which  is  used  to  generate  the  PRS:  K  -  2  -  1;  hence, 

a  typical  PRS  might  have  a  period  of  3,  7,  15,...,  1023,  2047,  4095,  etc. 
Figure  2.4  illustrates  two  typical  pseudorandom  sequences  of  period 
15  and  31,  along  with  their  respective  autocorrelation  sequences.  Note 
that,  for  i  +  nK,  for  integral  n,  *  -1. 


16 


24 


32 


40 


48 


56 


(a)  K-15 


Figure  2.4  Typical  PRS  and  autocorrelation 
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Construction  of  the  FLAC  transmitter  is  shown  schematically  in 
Figure  2.5. 


Figure  2.5  Discrete-time  FLAC  transmitter 
The  output  is  seen  to  be 

m  m  k  ^  *  Oj  if  2, . . . 

where  the  sequences  {bfc}  and  {m^}  are  the  original  information 
bearing  sequence  and  the  resultant  encoded  sequence,  respectively,  both 
of  which  were  introduced  earlier. 

At  the  receiver,  the  analog  input  is  sampled  every  T  seconds  to 
produce  an  input  sequence  {R^}.  Thus,  under  Dopplerless  conditions 

each  transmitted  bit  m  ,  s,  is  sampled  N  times  (N  will  usually  be 

[-1  k 
lKJ 

taken  as  an  integer),  with  each  sample  shifted  into  a  shift  register 

train  (Figure  2.6)  composed  of  two  individual  shift  registers,  each 

.th 


of  length  L.  At  the  P 


aP  =  aP(nT)  *  R  R 
n  n-Pn-L-P 

Y  is  taken  to  be 
n 


bit  position  of  each  register  the  product 

is  calculated,  and  the  receiver  output 
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l-i  p  L;L 

Y  =  y  a  -  y  r  ..  r  T  p 

n  poo  n  p-o  n  n" 


(2.21) 


In  order  to  obtain  an  expression  for  Y  in  terms  of  s,  note  that, 
with  Doppler  C,  and  under  noiseless  conditions,  (recall  that  an  =  a(nNT)), 


R  =  R(nT)  -  a(£nT)  -  m  s(CnT) 

n 


(2.22) 


m 


Since  K  >  1  the  whole  part  of  5n/N  may  be  used  in  the  subscript  of 

m,  and,  for  the  signal  being  considered,  s(5nT)  -  s  ,  so  that  (2.22) 

[  N  1 

may  be  written  as 

„  o  (2.23) 

R  =  m  s 

n  riEi  riH] 
lNKJ  1 N  J 

By  defining  and  sk  as  zero  for  k  <  0,  equation  (2.21)  may 

be  written  as 


L-l 


\  "  J  m.E(n-P),  mr£(n-P)-q1  B(5(n -  P)T)s(5(n -  P  -  L)T) 
P"°  1  KN  J  1  KN 

(2.24) 
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for  all  n  >_  0.  Equation  (2.24)  is  the  general  form  for  the  receiver 
output;  for  the  PRS-type  signals  being  used  (2.24)  becomes 
L-l 

Yn  %L  °trS(n-P)1  mfS(n-P)-^L  S,5(n-P)  S  g(n-P)-SL. 

1  KN  J  1  KN  JlNJl  N  J 

(2.25) 

If  the  m's  in  (2.25)  are  equal  for  0  <_  P  <_  L  -  1,  it  is  desirable  to 
maximize  Y^  for  5  =  It  in  which  case  the  shift  register  length  L 
should  be  chosen  as  L  =  KN.  For  if  L  is  chosen  in  this  way,  the 

subscripts  of  s  in  (2.25),  for  5  =  1,  become  [-■  -— ]  and 

N 

n  -  P  N  -  P 

l — « - K]  ■  [ — — — ]  -  K,  thus  assuming  a  peak,  as  can  be  seen  from  (2.20). 

N  N 

A  typical  output  sequence  is  illustrated  in  Figure  2.7  for  K  *  15, 

N  =  5,  L  -  KN  =  75,  and  5  =  1,  with  the  information  sequence  b^  as 

Indicated  thereon.  A  feeling  for  the  effects  of  Doppler  (for  the  same 

values  of  K  and  L)  can  be  obtained  from  Figure  2.8,  a  graph  of 

N  Y2L  (  the  2L  was  chosen  arbitrarily)  versus  5,  with  m^  assumed 

to  be  +1  for  all  k.  Although  any  plot  of  Y  versus  5  is  very 

n 

dependent  upon  the  PRS  employed,  and  upon  the  shift  register  length  L, 
some  generalizations  can  be  gleaned  from  Figure  2  8.  First,  note  that 
the  value  of  Y^^  becomes  very  erratic  and  unpredictable  for  5  out¬ 
side  the  range  of  approximately  (K  -  1)/K  to  (K  +  1)/K,  and  secondly, 
the  value  of  N  does  not  seem  to  have  a  drastic  effect  upon  the  values 
in  the  graph.  Specifically,  the  larger  value  of  N  seems  to  offer  no 
improvement  in  system  performance.  In  fact,  suppose  N  is  taken  to 
the  limit  in 

l  l  NK"1 

Y '  =  -jr  Y(nT)  =±  £  s(5(n  -  P)T)s(5(n  -  P  -  L)T) 

n  ”  ”  P-0 

For  convenience,  the  product  NT  may  be  assumed  constant  at  NT  =  1, 
corresponding  to  the  transmission  of  one  bit  of  the  PRS  per  second,  so 
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that,  letting  t  =  nT,  there  obtains  [7  ] 

Y'(t)-lim  ^  y(nT)  -  /  sU(t  -  T))s(S(t  -  t)  -  SOdx  (2.26) 

N  -*• "  0 

which  is  the  counterpart  of  equation  (2.8),  with  =  +1  for  all  k. 
If  (K-l)/K,  or  5  -  (K  +  l)/K,  the  value  of  the  integral  in  (2.26) 

is  approximately  zero,  which  indicates  that  the  system  is  usable  only 
for  1  -  1/K  •-<  F  «  1  +  1/K,  a  region  which  shrinks  with  increasing 

K. 

The  upshot  of  the  above  argument  is  that  performance  of  the  FLAC 
receiver  with  regard  to  Doppler  cannot  be  dependably  improved  by 
sampling  at  a  higher  rate  at  the  receiver,  a  result  which  is  to  be 
expected  from  the  analysis  undertaken  previously. 
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C.  Interference  Sensitivity:  The  FLAC  System 
Suppose  there  are  two  signals  impinging  upon  the  FLAC  receiver, 
one  of  which  it  is  desired  to  detect,  and  the  other  of  which  is  con¬ 
sidered  to  be  interference.  Denote  by  Sd»  K^,  and  Nd  the  parameters 
of  the  desired  signal,  and  by  and  ^  the  parameters  of  the 

undesired  signal,  where  these  symbols  play  the  same  role  as  their 
unsubscripted  counterparts  in  Section  II-B.  Further,  let  the  desired 
signal,  as  transmitted,  be  given  by 


a,  *  m  k  s  -  m  .  s(t) 

ir1  [k  1 

Kd  a 


kN^T  <_  t  <  (k  +  DNdT»  k 
be  given  by 


for 


and  let  the  undesired  signal 


6k  J  k  , 
k  [-^-1 

lK  J 
u 


X  *  J  ,  X(t)  for 
k  [_£_] 

V 


kNT<t<  (k  +  1)  N  T.  Then  the  received  sequence  OO  is  given  by 
u  —  u 


R  *  m 


r  V 

i8ridii 

+  J  _  A 

[ial]  [! 

VI 

9 

In  k 

N, 

N  K  J  l 

N  J 

u 

1  d  d 

d 

u  u 

ference  (2.25)  becomes 
L-l 


n 


l  mn 


P-0 


L-l 

+  l  mn 

P-0 


m 

5d(n-P-L)“ 

s 

s 

5d(n-P-L) 

KdNd 

KdBd 

Nd 

Nd 

. 

- 

Rd<-  - « I 

J 

\  (n-P-L)' 
u 

s 

fCd(n-p)' 

X 

\  (n-L-P) 
u 

KdNd 

K  N 
u  u 

ss° 

_ i 

N 

u 

m  m 

L-l 

+  l  m 

P-0 


[Cd(n-P-L) 


KdNd 


L-l 

+  l  Jn 

P-0 


Jr 


C^n-P) 


K  N 
u  u 


Sftd(n-P-L) 


Xr, 


N. 


C^n-P) 


N 


Cu(u-p) 


K  N 
u  u 


Jr, 


(^(n-P-L) 


K  K 
u  u 


xr 


?u(n-P) 


Xr, 


C  (n-P-L) 
u  _ 


N 


(2.27) 
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The  first  and  fourth  terms  in  (2.27)  are  autocorrelation  sums,  whereas 
the  second  and  fourth  terms  are  crosscorrelations  of  the  two  signals  s 

and  X.  It  is  a  reasonable  assumption  that  the  periods  and  of 

s  and  X  are  unequal,  and  that  the  crosscorrelation  of  the  two  sig¬ 

nals  will  be  approximately  zero.  Therefore,  the  second  and  third  terms 
in  (2.27)  will  hereafter  be  neglected,  and  the  only  term  which  will 
constitute  interference  will  be  the  autocorrelation  of  x. 


Since  K  and  K,  are  of 
u  d 

BL 

the  form  2  -  1,  for  integral 

BL, 

the  relation  between  K  and 

u 

K,  must  be  either 
d 

K .  -  £  K  +  £  -  1, 
d  u 

£  =  2,  3,  ... 

(2.28) 

or  K  -  £  +  1 

V  rn  _ y - 

t  * 

£  “  2,  3,  ... 

(2.29) 

Supposing  this  relation  is  given  by  (2.28),  end  that  L  =  K^Nj,  a 
contribution  from  the  fourth  term  in  (2.27)  requires  that  the  subscripts 


of  X  differ  by  kKu>  k  *  0,  +1,  +2,  ....  which  is  satisfied  for 


k  N  K 
u  u 


kN 


’u  Nj (IK  +  £  -1)  *  £n“  ’  k  “  °’  2 ’ 


(2.30) 


where  the  approximation  holds  whenever  £K  >>  £  -  1.  Thus,  if  N^,  N  , 

N  U 

and  £  are  chosen  such  that  -rjf-  m  2,  for  example,  then  (2.30)  becomes 

£Nd 

£  »  0,  2,  4,  ....  which  is  impractical,  so  that  appreciable  interference 

u 

will  not  be  observed. 

If  K,  and  K  are  related  by  (2.29),  the  autocorrelation  of  the 
d  u 


interference  signal  contributes  to  the  output  for 

k£N  N 


k£N  K 

r  .  - - 

Sj  N,  (K  -  £  +  1) 


N 


£N 


d  •  u 

,2 


N, 


If  is  chosen  as  “2,  then 

£Nd  fcNd 


— -  •  2lZ  >_  8,  and  once  again  the  values  of  £  required  for  appreciable 
Nd 

interference  are  impractical. 


'  : 

From  the  above  discussion  it  is  apparent  that,  if  the  bit  rates 
NT  and  the  periods  K  of  all  transmitters  operating  within  a  con¬ 
tiguous  region  can  be  controlled,  then  the  FLAC  receiver  can  be 
"tuned"  to  a  particular  transmitter  by  properly  selecting  the  receiver 
shift  register  length  L,  and  the  sampling  period  T. 
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III  The  Alternating  Forward/Reverse  Sequence  Autocorrelation  Receiver 
A.  Continuous  Time 

In  this  section,  the  Alternating  Reverse  Sequence  Autocorrelation 
(ARSAC)  receiver  is  discussed.  The  ARSAC  system,  conceived  by 
Cornett  [  8],  assumes  the  transmitted  signal,  x(t),  is  given  by 


x(t) 


m  t 

[h 


s(t). 


(3.1) 


where  t 


(t)2T  +  2 


^fc)2T 


(T  -  (t)2T),  and  is  a  differen¬ 


tially  encoded  version  of  some  information  sequence  {b^}^^  as  *n 
(2.6).  The  functions  £t  and  [^jp]  are  illustrated  in  Figure  3.1 
for  £  e  (.75,  1,  1.25}.  Note  that  the  use  of  t  as  the  argument 
for  s ( • )  in  the  transmitter  has  the  effect  of  transmitting  s(.) 
alternately  forwards  and  backwards.  Denoting  the  received  process  by 
r(t),  the  ARSAC  receiver  output,  y(t),  is  given  by 
T 


y(t) 


r (t  -  x)r(t  -  2T  +  i)di. 


(3.2) 


Suppose  that  the  signal  component,  S(t),  of  the  received  process, 
r(t)  ■  S(t)  +  n(t),  is  given  by  S(t)  *  x(t),  and  that  n(t)  is  GWN. 
Defining  z(t)  *  E{y(t)},  we  obtain,  for  t  ■  kT: 


z(t=kT) 


\-l  \  -  2 


T 

» 

s((kT-  T))s((kT-  2T  +  i))dT 

0 


“k-l  2 


s((k-u)T)s((k-2  +  u)T)du 


(3.3) 


0 

Suppose  that  there  exists  an  a,  b  such  that  for  all  a  e  [a,  b)  we 
have  [a]  ■  h,  h  a  constant.  Then  if  h  is  even,  we  have  aT  *  (a-h)T 
and,  if  h  is  odd,  we  have  aT  ■  (  1  +  h  -  a)T.  Consequently,  we  find 
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from  (3.3)  that  for  k  even,  (k  -  u)T  ■  (k  -  2  +  u)T  =  uT  and  for 
k  odd,  (k  -  u)T  =  (k  -  2  +  u)T  -  (1  -  u)T,  for  all  u  e  [0,  T).  Finally, 
we  find  that  (3.3)  may  be  written  as 


z(t  -  kT)  -  mk_  1  \_2  cs(°). 


(3.4) 


where 


CS(Y)  *  I  8(T)8(T  -  Y)dT. 


We  note  that  (3.4)  and  (2.4)  are  identical,  indicating  that  the  operation 
of  the  FLAC  system  and  the  ARSAC  system  are  identical  at  t  *  kT, 
when  there  is  no  Doppler  shift. 

Now,  suppose  that  the  signal  component  of  the  received  process  is 
given  by 


S(t)  -  m  s  (Ct). 

i£i  d 


(3.5) 


Using  (3.2),  with  r(t)  -  S(t)  +  n(t),  n(t)  GWN,  and  z^(t)  ■  E{y(t)}, 
we  obtain 


i 

ZC(t>  *  \  ".Ut-2T+,>.  8d(5(t  -  T)).d«(t-2T  +  ,»dt. 

0  ^  t  T 


(3.6) 


Making  the  change  of  variable  u  =  -^jr,  and  defining  x  *  and 
z^(x)  ■  z^(t),  we  obtain 


z'  (x) 


T  f*  _ 

K  J  "tx-u]  m[x-2^u]  Sd((x-U 


)T)  sd((x- 2C  +  u)T)du. 


(3.7) 


Partitioning  the  interval  [0,  C]  as 


0  “  a0  <  al  <  a2  <  *  '  ,<:  aL  ”  ^ 
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such  that,  with  x  =  k  +  (x)^,  we  have  [x  -  u]  ■  k  +  n^,  and 

[x  -  2?  +  u]  =  k  +  P^  for  all  u  e  [a^^,  a^),  we  see  that  (3.7)  may 

be  written  as 


z'(x) 


where 


r  £  mk  +  n  mk  +  P 
4  £"1  £  K  *£ 


sd((x-u)T)sd((x- 2?  +  u)T)du,  (3.8) 


£  -  1 


((x)^  -  u  -  n^)T,  if  k  +  n^  is  even 


(x  -  u)T  ■< 


(1  -  (x)^  +  n^ +  u)T,  if  k  +  n^  is  odd; 


(3.9a) 


and 


(x-25 


J((x)x  -  25  +  u-P^T,  if  k  +  P^  is 


even 


(3.9b) 


(1-  (x). +2£  +  P  -u)T,  if  k  +  P  is  odd 

V.  1  *  Jt 


Table  3.1  gives  all  the  parameters  needed  for  the  computation  of  (3.8) 
and  (3.9)  for  1/2  <  5  <  3/2.  The  maxima  and  minima  of  (x)  should 


occur  for  those  values  of  x  for  which  (x  -  u)T  -  (x-2£  +  u)T  in  (3.8). 
After  investigating  the  possibilities  with  the  aid  of  (3.9),  it  is 
found  that  the  maxima  and  minima  will  occur  when 


5,  1/2  <  5  <  1 


(x)x  -  *<£  -  1,  1  <  £  <  3/2 


0,  5  -  1 


Using  Table  3.1,  for  1/2  <  £  <  1  and  (x).  ■  £,  we  find 

ft 


z£(x)  ml\\- 1  i 


J  s , ((£  -  u)T)du,  for  k  even 
0  d 


/  s,((l-  £  +  u)T)du,  for  k  odd 

U)  d 


(3.10) 


Table  3.1 


Parameter  values  for  the  computation  of  the  output  of  the  ARSAC  receiver 

(a)  1/2  <  £  <  1 


(x)1  e  {.} 

£ 

ft 

a£  -  1 

a£ 

[0,  5  -  1/2) 

1 

0 

-2 

0 

(x)x 

2 

-1 

-2 

(x)^^ 

25  -  M1  -  1 

3 

-1 

-1 

25  -  (x)x  -  1 

5 

[5  -  1/2,  2K  -  1) 

1 

0 

-2 

0 

25  -  (x)x  -  1 

2 

0 

-1 

25  -  (x)x  -  1 

(x^ 

3 

-1 

-1 

(x^ 

5 

[25  -  1,  5) 

1 

0 

-1 

0 

(x^ 

2 

-1 

-1 

(x^ 

5 

[5,  1) 

1 

0 

-1 

0 

25  -  M1 

2 

0 

0 

25  -  (x)x 

5 

(x).  e  {.} 

£ 

(b).JL 

-  1 

a. 

a„ 

[0,  1/2) 


[1/2,  1) 


[i, 
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Similarly,  for  1  <  5  <3/2  and  (x)  -  f.  -  1,  we  obtain 

fc-i 


s,  ((£  -  1  -  u)T)du,  for  k  even 
a 


zc  (x)'?mk\-3  V 


s,  ((2  -  £  +  u)T)du,  for  k  odd 
a 


0.11) 


s,  ((1  -  C  +  u)T)du,  for  k  even 
a 


s  ((£-u)T)du,  for  k  odd 
a 


Now,  defining  the  gain  of  the  system  against  Doppler,  G,  as 


G  *  |z' (x) |/C  (0), 

6  d 


P  =  T  •  max  s ,  (uT) , 
max  rA  1Nd 

ue[0,l) 

we  obtain,  for  1/2  <  £  <  1.  and  (x)^  *  £» 


G  >  i  +  _Eax_  • 

G  -  £  £  C  (0)  ’ 

sd 


(3.12) 


whereas,  for  1  <  £  <  3/2,  (x)^  =  5-1,  and  m^_  3  “-m^  ^  mjt_  2’ 


r  1  .  (1  ~  C)  max 
°  -  C  +  5  C„  <0> 

Sd 


(3.13) 


The  lower  bounds  of  (3.12)  and  (3.13)  are  illustrated  in  Figure  3.2  for 

1/2  <  5  <  3/2  and  P  /C  (0)  e  {1.5,  2}. 

max  Sj 


B.  The  Discrete  Time  ARSAC  Receiver 

As  before,  the  carrier  signal  s(t)  will  be  assumed  to  be  a 
sampled  and  held  version  of  a  PRS  of  period  K,  but  with  this  system, 
the  first  K  bits  of  the  sequence  must  be  transmitted  alternately 


forwards  and  backwards.  By  defining  t„  M  as  t  =  J  +  [-]  (2K  -  2J  -  1) , 

tv  j  W  K  j  N  K 

with  J  *  (■£■]  modulo  2K,  the  carrier  signal  may  be  taken  as 
N 


s  ( t)  -  s 


(t/T) 


which  will  be  shortened  to 


when  no  confusion 


will  result.  Here,  N  and  T  have  the  same  meaning  as  in  the  previous 
section,  so  that  one  bit  of  the  PRS  is  transmitted  every  NT  seconds. 

A  graph  of  the  function  £ t  versus  t  can  be  seen  in  Figure  3.3,  for 
5  e  {.75,  1,  1.25}. 

The  ARSAC  transmitter  is  illustrated  schematically  in  Figure  3.4. 

The  transmitter  output  is  ^t/T*  w^ere  denotes  a(t)  for 

kNT  t  <  (k  +  1)NT,  k  =  0,  1,  2...,  and  tR  denotes  [^] . 

At  the  receiver,  the  input  is  sampled  every  T  seconds,  providing 
an  input  sequence  {R^},  with  each  sample  entered  into  a  shift  register 
train  composed  of  two  shift  registers,  each  of  length  L,  one  of  which 
shifts  to  the  right,  and  the  other  to  the  left,  as  is  illustrated  in 
Figure  3.5. 

The  receiver  output  Y  is 

n 

L-i  p  1—1 

Yn  % l0  Sn  =?lQ  Rn  -  P  Rn  -  2L  +  P  +  1 

With  Doppler  5,  R  =  «(CnT)  =  a  .  =  m,_  *  S7— ,  so  that  Y  is  given  by 
*  n  (Cn) ^  £n’  n 


Yn  m(C(n-P))v  m(C(n-  2L  +  P  +  1))  S£(n-P)  SC(n-2L  +  P  +  l) 


(3.14) 


jo 


Equation  (3.14)  is  valid  for  all  n>0  if  it  is  agreed  that  m^-s^O 
for  k  <  0. 

Figure  3.6  illustrates  a  typical  output  sequence  for  K  =  15,  N  =  5, 

L  -  KN  -  75,  and  K  =  0.81.  Note  that  the  outrut  becomes  decorrelated 
between  the  peaks,  which  are  manifested  distinctly  above  what  might  be 
termed  correlation  "noise".  This  represents  another  advantage  over  the 
FLAC  system,  in  that  estimates  of  the  information  bits  may  be  made 
directly  from  the  output  amplitude,  whereas  such  estimates  must  be 
made  at  specific  times  when  using  the  FLAC  receiver. 

The  positions  of  the  peaks  in  time,  and  the  relative  spacing 
between  them  may  be  determined  somewhat  houristically  by  the  following 
argument:  After  operation  of  the  receiver  has  begun,  the  first  data 
sample  reaches  the  (L  -  l)st  position  of  the  first  shift  register 
after  L  -  1  steps.  Thereafter,  an  entire  PRS  group  of  K  bits 

passes  through  this  position  and  into  the  second  shift  register  approxi¬ 

mately  every  NK/£  steps.  Thus,  one  might  expect  that  a  peak  will 
occur  in  the  output  sequence  Yn  whenever  n  is  an  integer  close  to 
L  -  1  +  kNK/5,  k  -  1,  2,  3,....  It  will  be  shown  below  that  this  integer 
is  actually  n  -  {L  -  1  +  where  {x}  denotes  the  smallest  integer  ^x. 

In  order  to  visualize  the  effects  of  Doppler,  assume  for  the 
moment  that  =  +1  for  all  k  and  define  An  as 

An  ■  L  —  1  +  ^  _  n,  k  1,  2,  3,..., 

where  n  is  taken  as  n  «  (L  -  1  +  kNK/U-  Then  -1  <  An  <  0,  and  the 

subscripts  of  s  in  (3.14)  become 

g(n  -~pj  =  5(L  -  1  +  An  -  P)  -  kNK  +  £(L  -  P)  -  5(1  +  An) 


and 


i+«  q 


Figure  3.6  Typical  ARSAC  Output 
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Figure  3.3  indicates  that 

kNK  -  x  =  kNK  +  x 
for  all  x  such  that 

mN  <x<  (m+l)N,  m  “  0,  1,  2,..., 

(equality  may  hold  on  the  left  for  m  =  0) ,  so  that  the  subscripts  in 
(3.14)  will  be  equal  for  P  such  that 

mN  <  £(L  -  P)  -  5(1  +  An)  <  (m  +  1)N 

and 

mN  <  £(L  -  P)  +  5An  <  (m  +  1)N. 

Choosing  the  tighter  upper  and  lower  bounds  for  5(L  -  P)  begets 
mN  +  C(1  +  An)  <  5(L-  P)  <  (m  +  l)N  -  5An,  -1/2  <  An  <0 

and 

mN-?An  <  £(L  -  P)  <  (m  +  1)N  +  5 (1  +  An)  ,  -1  <  An  <  -1/2 


whence 

L  +  An-  (m  +  l)|<  P<L-l-An-n>|  ,  -l/2<  An  ^0  (3.15) 

L-  1  -  An-  (m  +  l)7<  P  <  L+  An-  -1<  An  1-1/2  (3.16) 

The  rationale  for  the  particular  choice  of  n  may  be  interjected  here. 
For  consideration  of  equations  (3.15)  and  (3.16)  indicates  that  the 
greatest  possibility  oi  satisfying  these  inequalities  occurs  whenever 
An  *  -1/2.  Since  the  difference  between  the  maximum  and  minimum  values 
of  An  must  be  at  least  one,  the  optimum  range  of  An  is  -l<An<0, 
as  was  chosen. 

To  insure  the  existence  of  an  integer  between  the  bounds  in  (3.15) 
and  (3.16)  for  every  m  requires  that  the  segment  width  defined  by 
these  bounds  be  greater  than  1,  which  in  turn  requires  that  N  >  25- 
Thus,  if  the  maximum  possible  value  of  5  is  denoted  as  Cmax.  it  is 
necessary  that  N  >  25max  for  reasonably  good  operation  of  the  system. 
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Note  that  the  larger  the  value  of  N,  the  more  integers  will  lie 
inside  the  segments,  with  resulting  better  system  performance. 

The  above  discussion  indicates  that  bit  misalignments  will  occur 
at  the  Pth  position  of  the  correlator  shift  registers  whenever  P 
satisfies 

L  -  1  -  An  -  —  <  P  <  L  +  An  -  1117,  -1/2  <  An  <_  0  (3.17) 

E,  —  —  s 

or 

L  +  An-rn^-  ^  P  1  L  -  1  -  An  -  ,  -1  <  An  £  -1/2  (3 .18) 

tor  „  -  1,  2,  3 .  [f  a  +*■>)],  -1/2  The  lengths  of  the 

intervals  in  (3.17)  and  (3.18)  are  given  by 
0  <  1  +  2 An  <_  1,  -1/2  <  An  £  0 

and 

0^-1  -2An  <1,  -1  <  An  1  -1/2 

Hence,  at  most  two  integers  may  lie  within  the  interval  defined  by  (3.17) 
so  that,  for  the  worst  possible  case,  a  total  of 
2[|(L  +  An)]  (for  -  \  1  An  5.0) 

bit  misalignments  may  exist  in  the  shift  registers. 

Peak  degradation  will  occur  due  to  these  misalignments  only  if  the 
appropriate  adjacent  bits  of  s  are  of  opposite  polarity,  as  is  shown 
in  Figure  3.7,  an  illustration  of  a  small  section  of  the  correlator 
shift  registers. 


No  Degradation^  Degradation  Occurs 

Figure  3.7  Bit  misalignment  causing  peak  degradation 


If  degradation  does  occur,  the  magnitude  of  the  peak  will  be  lessened  by 
2  (+1  is  subtracted  instead  of  added)  for  each  bit  which  is  misaligned 
and  of  polarity  opposite  that  of  an  adjacent  bit.  Letting  A^(R) 
denote  the  maximum  number  of  polarity  reversals  in  R  bits  of  a  PRS 
of  period  K  (with  successive  groups  of  K  bits  transmitted  alternately 
forwards  and  backwards),  the  maximum  amount  of  degradation  AY^  which 
may  occur  due  to  Doppler  is  (2  bits  misaligned  per  reversal)  x  (2  sub¬ 
tracted  per  misaligned  bit)  x  (number  reversals),  so  that 

AY.  =  4A,,(UL/N})  (3.19) 

d  K. 

The  value  {£L/N}  is  simply  the  number  of  PRS  bits  which  ' fit"  into 
each  of  the  shift  registers  whenever  a  peak  occurs. 

As  an  example,  suppose  K  =  7,  and  that  the  PRS  sequence  is  trans¬ 
mitted  as 


r-l,  1,  -1,  -1,  1,  1,  O  / 1 ,  1,  1, 


S6  S5  ‘ 


Supposing  further  that  L  =  NK,  and  K  -  0.8,  then  (3.19)  is  evaluated 
as 


AYd  -  4A?  ({0.8  •  7})  -  4A?  ({5.6})  -  4A?  (6). 

In  the  sequence  sQ,  s1>...,s5  there  are  3  polarity  reversals,  whereas 
in  the  s,,  sc,...s.  there  are  only  2  such  reversals.  Thus  A  (6)  -  3, 
and  (3.19)  evaluates  to  AY^  ■  12. 

In  order  to  prove  the  advantages  of  the  ARSAC  receiver  with  regard 

to  Doppler,  let  the  normalized  peak  amplitude  Y'  be  defined  as 
L  -  AY 

Y' - - 

N 

Now,  if  L  is  proportional  to  N,  i.e.,  L  ■  BN,  then  Y'  becomes 

AY.  4A  ({£0}) 


N 
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which  tends  to  B  as  N  tends  to  infinity,  independently  of  £.  Thus, 
relative  peak  degradation  due  to  Doppler  (with  =  +1  for  all  k, 
as  was  assumed  earlier)  can  be  nade  arbitrarily  small  in  the  ARSAC 
system  by  increasing  the  sampling  rate  at  the  receiver. 

With  the  restriction  that  =  +1  removed,  the  above  analysis 
may  have  shortcomings,  since  the  value  of  Y  ,  for  n  =  (L-l  +  kNK/C), 
may  depend  not  only  upon  ra^  and  m^  _  ^ ,  but  also  upon  mjc  +  ^  anc* 
m^  2-  The  continuous  time  analog  of  this  problem  was  treated  in 
Section  IV-B. 

An  analytical  description  of  the  situation  is  achieved  by  con¬ 
sidering  the  subscripts  of  m  in  (3.14),  the  first  of  which  becomes, 
for  n  as  given  above, 

(S(n-P))K=  (kNK+aL-P)  -  5(l  +  An))K 

Employing  the  definition  of  tR,  this  equation  becomes 

«<■  -  p»k-  [*"  +  P>  ~  ^1  -  *  *  lt<L~  *M1 

Hence,  the  subscript  under  consideration  remains  constant  at  k  if 

0  <  j(L-P)  -  ^l.+i"),  1 
-  KN 

N 

which  is  satisfied  for  L  -  1  -  An  -  K~  <  P  L  -1. 

Similar  treatment  for  the  second  subscript  of  m  indicates  that 
(5(n  -  2L  +  P  +  1))„  remains  constant  at  k  -  1  if 

i\ 

L  +  An  ••  k|  1  P  1  L  -  1 , 

so  that,  since  -1  <  An  £  0,  both  inequalities  are  satisfied  for 

{L  -  k|)  <_?  <_  L  -  1  <3  -20) 

For  P  outside  this  range  the  two  subscripts  become  k  +  1,  and  k  -  2, 
respectively.  Therefore,  under  the  worst  possible  conditions,  which 
correspond  to  the  products  ra^  m^  _  ^  and  m^  +  ^  “'k  -  2 


A 


* 
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opposite  polarity,  and  to  the  exact  alignments  of  all  PRS  bits  in  the 
Pth  position  of  the  shift  registers,  for  0  <  P  <  {L-Kj}  = j ,  the 
peak  amplitude  is  degraded  by 


AY 

m 


2<J  +  1), 


l 


j  >  0 
j  <  0 


Consideration  of  equation  (3.20)  shows  that,  if  L  is  such  that 

L  -  KN/(  <  0,  then  the  degradation  AY  will  be  zero.  Thus,  L  can 

m 

be  chosen  so  that  L  <  KN/£,  which  is  satisfied  for  L  =  [KN/£  1. 

max 

If  L  is  chosen  in  this  manner,  only  the  degradation  due  to  bit 
misalignments  will  contribute  to  peak  amplitude  reduction. 

In  the  above  analysis  it  has  been  assumed  that  the  received 
sequence  {R^}  is  composed  of  signals  emanating  from  a  single  trans¬ 
mitter,  and  that  no  interference  is  present  from  other  systems  in 
operation  near  the  receiver.  If  such  is  not  the  case,  system  operation 
may  be  impaired  to  the  point  that  the  desired  signal  is  totally 
undetectable.  Consideration  of  this  problem  and  of  methods  for  achieving 
its  remedy  will  be  undertaken  in  the  following  section. 
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C.  Interference  Sensitivity:  The  ARSAC  System 

For  convenience,  let  it  be  assumed  at  the  outset  that  the  encoded 
message  bits  are  +1  for  all  time  for  both  the  desired  and  the  undesired 
signal.  Then  the  transmitted  signals  are  given  by 


and  x; 


da  u  u 


and  (3.14)  becomes 


L-l 

Y  *  l  s 

n  A- 


“  pio  <5d<"-P>>KJ,N/(«dt''-2L  +  P  +  1)>K 

d  d 


xd  *  d 


L-l 

X<"  -  P))K,,N  X<"-2L  +  P  +  1»K  ,N 

d  d  u  u 

L-l 

"p-0  S(Cd(n  -  2L  +  P  +  1»  .  >N 

da  u  u 

L-l 

+P-0  X(”  -  P))K  ,N  X("  -  2L  .  ?  .  i)TK  N 

U  U  U  U 


(3.21) 


Again  it  will  be  assumed  that  the  second  and  third  sums  are  negligible, 

so  that  only  the  first  and  fourth  terms  contribute  to  Y  .  No  further 

n 

analysis  is  required;  from  the  results  of  Section III-B  one  sees  that 

the  fourth  term,  which  corresponds  to  the  interference,  peaks  at 

n  *  (L  -  1  +  k  K  N  /£  }, 
u  u  u 

for  k  m  1,  2,  3,....  Furthermore,  the  maximum  possible  value  of  this 

term,  L,  Is  the  same  as  that  for  the  first,  desired  term.  Hence,  the 

peak  due  to  the  desired  signal  and  the  peak  due  to  interference  may  be 

indistinguishable  in  amplitude,  and  since  and  £  are  not  in  gen- 

d  u 

eral  controllable,  the  relative  position  in  time  of  these  two  peaks 
cannot  be  dependably  established. 
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Although  the  above  results  were  deduced  with  the  assumption  that 
the  message  bits  for  each  signal  were  constant  at  +1,  an  unreasonable 
assumption,  it  is  apparent  that  an  equivalent  situation  may  occur  if 
a  string  of  +l’s  or  -l's  occurs  in  the  message  bit  sequence  corresponding 
to  the  interference  signal.  Therefore,  it  must  be  concluded  that  the 
ARSAC  system  exhibits  poor  interference  characteristics. 


45 


IV.  Noise  Performance 

In  this  section  we  investigate  the  noise  performance  for  the  family 
of  predistorted  replica  correlation  receivers.  We  will  be  concerned 
with  the  computation  of  the  cumulative  distribution  function  (cdf)  of 
the  statistic  z,  given  by 


1  _T  _ 

z  '  H  R1  S2- 


(4.1) 


where 


Ri  "  Si  +  Ni 


Ri  =  ^ri,l*  ri,2* 


...  r±>N). 


Si  "  (si,l»  Si,2’  **•’  Si,N)* 


Ni  =  (ni,!’  ni,2’ 


n  ,  i  c  (1,  21,  l  c  (1,  2,  ...,  N}  is  normally  distributed  with  mean 

2  2 
zero  and  variance  (we  will  use  the  notation  n^  ^  N(0,  o^)),  and 

E  {n  n  ,  }  =  o2  6,  .  6,  ,  .  S.  is  the  (deterministic)  signal  com- 

i,i  j,k  i  i,j  k,l  i 

ponent  of  the  received  data  vector,  R^.  Several  solutions  to  similar 
problems  are  available  in  the  engineering  literature  (e.g.  see  [9  ]  -  [11]); 
however,  since  the  results  treated  there  are  not  applicable  to  the  problem 


at  hand,  they  will  not  be  dealt  with  here. 

We  begin  our  analysis  by  computing  the  characteristic  function 


ij)  (5),  defined  by 
z 

<j>z(0  =  E  {eUz},  (4-2) 

where  E{ • }  denotes  statistical  expectation.  Feller  [12]  suggests  that 
if  and  X£  are  independent  random  variables  with  cdf's  and 

F2,  and  characteristic  functions  (ch.f.'s)  ^  and  $2*  respectively, 
then  y  =  x^x2  has  the  ch.f.  <t>^  given  by 
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ye) 


<t>2  (Cx)dF 


oo 

l(x)  -  j 


(Cx)dF2(x). 


(4.3) 


12  2 

Making  use  of  (4.3)  and  the  fact  that  4)  (O  =  exp{j  s  £  -  y  o  £  ) . 

ri,k  1,K  1 

we  find  that  the  ch.  f.  for  the  random  variable  ^  s  rj  jj  r2  ^ 


given  by 


4>  (5)  =  exp 


V  % 


r  .1,2  2  2  2  2 
s„  ke  -  y  (o2  Sl  k  +  °1  s2  k)C 


2  2  2 
l  +  o1  o2  £ 


} 


'  ^7  J  exp 

i  -» 


-(J+Ol  o2  S2)  I  s1>|t  + j  °2.k  ° 


20 


(«  -1J 


2  2 

1+°1  °2  5 


or 


. .22  .2, -1/2  P  Sl,k  2  (°2  Sl.k  +  °1  S2,lP  } 

<«  "  <  1+°1  °2  5  >  H ; - f~2^ - - 

;  i  1  +  Oj^  o2C  J 


Now,  defining 


(4.4) 


C1  =  S1S2’ 


2  T  2  T 

C2  =  °2  S1  S1  +  °1  S2S2  * 


(4.5) 


and  noting  that  our  independence  assumption  implies  that 
N 

<J>Nz(e)  =  11  <f>  (5).  we  have 

k*l  yk 


_  N 

<f>  (N  0  =  (1  +  0?  o?  C2)  2  exp 
z  l  i 


P  Cl*  -  2  C2g  1 

CXPl  i  2  r2  I 
L  1  -h^  o2  £  J 


(4.6) 


since  a,  (^)  =  <h  (£).  Defining  the  standardized  random  variable,  p,  by 
Nz  N 


z 

Z  -  u 


(4.7) 


where 


UZ  *  E  U1"  T' 


and 
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2  2  °i  tt\  +  c2 

Eiz  }  -  y  =  — n +  -7, 

2  N  N2 


we  obtain 


♦n(f)  =  (1  +C2)  2  ex 


3  1  2 

fJ  V  -  2  B2  5  1 


i  1  + 


(4.8) 


where  8  =  a  02/Noz,  S1  =  S/°l°2*  and  S2  =  C2^°ia2*  Unfortunately» 
the  ch.f.  (4.8)  is  not  easily  inverted  to  obtain  the  denisty  function; 

however,  all  of  the  information  about  the  density  function  is  contained 
in  (4.8).  From  (4.8)  we  can  obtain  all  of  the  moments,  y,  =  E(n  }, 


using  the  relationship 


,  VV0! 

u.  =  (-j)  Hr—  I 

k  dt  t=o 


(4.9) 


Computationally,  a  somewhat  loss  tedious  approach  is  to  first  compute 


the  cumuiants  K^,  using  the  relationship 


log  ♦  <t)  -  1  1^  iif-,  (4.10) 

k*l 

and  noting  that  the  moments,  y^»  and  cumuiants,  K^,  are  formally  related 


by  [13,  P.  318] 


el<p  *Sc  tT  1  '  J  “k  kT 


(4.11) 


Kendall  and  Stuart  [13,  p.  69]  have  tabulated  the  first  ten  moments  in  terms 
of  the  first  ten  cumuiants,  which,  for  =  0  and  ^  “  1  are 
repeated  in  (4.12) 


y2  *  1 


^4  *  K4  +  3 
y5  -  k5  +  iok3 


(4.12) 
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V,  =  K,  +  15K.  +  lOK^  +  15 
do  4  3 

U?  -  K?  +  21K5  +  35K4K3  +  105K3 

U  =  KQ  +  28K  +  56K.-K-  +  35K?  +  210K.  +  280K^  +  105 
88  6  53  4  4  3 

H9  -  K9  +  361^  +  84K6K3  +  126K5K4  +  378K5  +  1260K4K3 
+  280K3  +  1260K3 

p  -  K10  +  45Kg  +  120K?K3  +  210K6K4  +  630K6  +  126K2  +  2520K5K3 

+  1575K?  +  2100K.K!?  +  3150K.  +  6300K2  +  945  (4.12; 

4  4  3  4  3 

The  reader  will  recall  that  for  a  standardized  normal  random  variable, 
the  only  nonzero  cumulant  is  *  1  and  the  nonzero  moments  are  given 
by 

u'  =  .L?A).r-  £-123  (4  13) 


Knowledge  of  the  first  few  moments  of  n  will  enable  us  to  compute 
an  Edgeworth  expansion  [14]  for  the  probability  density  function  f^, 
from  which  an  approximation  to  the  cdf  is  easily  obtained.  We 

now  proceed  in  developing  an  expression  for  the  cumulants  K^.  Defining 
h(0  -  log  <*>  (C/B) ,  we  have  o  .  ~ 

n  N  je.r+^r 

h(C)  --ylog(l+jO  -f  log(l  -  jO - - - \ - .  (4.14) 

Z  1  +  5 

Using  standard  partial  fraction  techniques,  and  noting  that 


00  /  1  ,k+l 

l°g(l  +jO  “  l  ^ - , 

k-1  K 


we  obtain 


00 

h(o  -  l  [(jo2k  +  +  Bi^^2k+1J- 

k=l 


It  follows  that  the  cumulants,  K^»  are  given  by 


^k!  (^  +  |b2  ),  k  =  2,  4,  6,... 


B13lCk! , 


k  =  3,  5,  7, . . . 


with  Kj-0- 


(4.15) 


(4.16) 


(4.17) 
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Substituting  (4.17)  into  (4.12)  we  obtain 


yx  =  o 


w2  -  l 


W3  =  3^3 

A  M  1 

u4  =  ^4  +  262)  +  3 

5  S'  3 
u5  -  5^3  +  y  Bj.3 


^  =  6:e6(f  +  I62+I6l>+fe4(!  +  i62)+15 

7  Ml  7 '  S  7 '  ^ 

M?  =  7!61B  (1  +  I  +  y62)  +  + 

U8  =  8:e8[(i  +  f)(f  +  yB2)  +  K+K] 

+  y66(f  +  yS2  +  4*i>  +  ^4(f  +  ie2)  +  105 

p„  *  9!M9(i  +  ^rN  +  80  +  |s*)  +  ■^e137(i  +  ^  +  ^32) 


12 

9'  5  9'  1 

+  TV  +lt«i8 


p10  •  10:b1°I<!  +  K><!  +  5V#  +  I5  +  I«2  +  K1 


(A.  18) 


10! 


•N+ifi^.N  +2ft2  10^6  N  K  12 

4  +  IV  +  4  +  62  +  2B1]  +  ^  (6  +  2^2  +  IV 


+  ( “i  +  ig  )  +  945 

48  l4  ?2'  * 


The  Edgeworth  expansion  is  a  series  representation  of  a  standardized 
probability  density  function  in  terras  of  Tchebychef f-Hermite  polynomials. 


the  first  term  of  which  is  the  standard  normal  density: 

1  2 


«(x) 


1  *  2  x 

vll  e 


(4.19) 


The  motivation  for  considering  an  Edgeworth  expansion  rather  than  the 
formally  equivalent  Gram-Charlier  Type  A  series  lies  in  practical 
applications  for  which  only  a  small  number  of  terms  are  used  [14].  In 


standardized  measure,  the  first  few  terms  of  the  Edgeworth  expansion 
are  found  from  [14,  p.  229]  to  be 


f<*>  =  ♦<*>  -  [j7  K3  *«>]  +  *4  * 


(4>  +10k2*<6)  , 

(x)  +  6!  K3  *  (x)J 


-  ‘iff  V<5)(x) +  if  W(u) +  W  k3  *(9)(x)]  «-20> 


An  added  benefit  in  using  expansions  such  as  (4.20)  is  the  ease  in 

obtaining  an  approximation  to  the  cdf,  F(x),  in  fact,  with 
x 


-  /  *(x')dx'  , 

-<*>  X 


F(x) 


j  f(x’)dx’ 
—00 


(2) 

(x) 


1 


+  t-Lr  J  3)  +  10  2  (5)  , 
+  [4T4  *  (x)  6!  K3  *  (x)] 


r_L  v  ,,<*)  .  35  „  „  (6)  .280  3.(8)  . 

‘  5!  K5  *  (x)  71  K3K4  *  (x)  9!  K3  *  (x)]> 


(A.  21) 


or,  after  "simplification," 


F‘*>  *  '•'(x)  +  WV  -  36is5  +  l»/<?  +  -  if6!®6' 


+  x(384(|  +  fy2)  -  yf^B6) 


+  x4(-81B3  +  6B185  -  A581B7(|  +  ^2)  +  70B3B9) 

+  x3(84(|  +  p2)  +  5B*B6) 

4  5  7  N  1  3  9 

+  xh(-B1B  +  15B16  (^  +  ^62)  -  35BJB’) 

+  x5e^B6)  +  x6(-B187(|  +  fy2)  +  y^B9) 

+  x8(-|  B3B9)]. 


(4.22) 


At  this  point,  it  seems  worthwhile  to  interpret  the  parameters  we 
have  been  using.  We  will  assume  that  the  signal  energy  is  the  same  in 
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T  'i' 

each  channel  of  the  correlator,  or  -  •s-)^>2"  Interpreting  the  Input 

1  T 

signal  energy,  Eg  ,  as  jjS  S]t  we  find  that  8  is  given  by: 
in  1  * 


8  = 


N1*  /I  +  y  SNRi  * 


(4.23) 


2.  2 


where  we  have  defined  y  =  Yj^  =  (1  +  y)/y,  and  the  input  signal- 

2 

to-noise  ratio  as  SNR  =  E  /o  .  We  now  define  a  correlation  factor, 

1  s ,  1 

in 

CORR,  by 


CORK  = 


T 

flf2 

T 

S1S1 


(4.24) 


which  is  easily  shown  to  satisfy: 


-1  <  CORR  <  1. 


Finally,  we  obtain 

N  •  CORR  *  SNR, 


0, 


and 


(4.25) 


82  -  Ny1SNRi. 


A  suitable  performance  index  for  binary  antipodal  signaling  is 

2|u  |  2 | CORR |  •  N1/2  •  SNR 

P  (y»  CORR,  SNR  ,N)  =  =  -  —  1  •  (4.26) 

z  /y+  (  1  +  y)  SNRi 

The  above  performance  index  (4.26)  is  illustrated  in  Figure  4.1.  For 


a  fixed  reference  correlation  receiver,  with  y  =  0,  and  reference 

T  T 

signal  S  c,  such  that  S  ,S  .  =  S,S,  and 
ref  ref  ref  1  1 


a  = 


T 

S.S  . 
1  ref 

T  ’ 

S1S2 


(4.27) 


our  performance  index  becomes 
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p  (0,  CORR,  SNR..N)  =  2a|CORRj/N  •  SNR. . 
a  11 


(4.28) 


We  can  now  define  an  efficiency  factor,  EF,  as 


EF (y ,  a,  SNR  ) 


p(Y,  CORR,  SNRi,N) 


/SNR, 


i  P  (0,  CORR,  SNR.,  N)  — — — 

a'  ’  ’  l  a/y  +  (1  +  y)SNRi 

(4.29) 

The  above  efficiency  factor  (4.29)  is  shown  in  Figure  4.2.  EF(y,  a,  SNR^) 
provides  us  with  a  strong  indication  of  the  tradeoffs  involved  in 
choosing  between  a  predistorted  replica  correlation  receiver  and  a  fixed- 
reference  correlation  receiver. 

With  the  above  interpretations  in  mind,  we  return  to  the  Edgeworth 
approximation  of  the  cdf  given  in  (4.22).  Since  H'(x)  is  the  well- 
tabulated  standard  normal  integral,  we  need  only  concern  ourselves  with 
the  non-normal  component  of  the  cdf,  H(x),  defined  by 


H(x)  -  F(x)  -  x|/(x). 


(4.30) 


The  Edgeworth  approximation  of  (4.30)  given  by  (4.22)  is  shown  in 
Figure  4.3  for  |x|  ^  5  and  several  choices  of  N,  SNR^,  y,  and  CORR. 
Only  positive  values  of  CORR  were  used  since 


H(x)  |  =  -H(-x) | 

-CORR  CORR 


(4.31) 


Figure  4.3  suggests  that  H(x),  the  non-normal  component  of  the  cdf  F(x) 
is  not  of  large  magnitude,  and  tends  to  decrease  with  increasing  N. 
Consequently,  an  investigation  of  the  behavior  of  $  (£)  ^or  lar8e  N 


would  seem  appropriate. 
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From  (1.8),  we  have 


where 


.(C)  =  (1  +  PT)  exp 


2  2-f  -JS1b3^  ■  2 

C  )  exp  < - 

1  +  6  € 


fi2  _  1  ....  _ *  0 

0  N(1  +  fjSNR  )  N  -*■  00  ’ 

CORR  •  SNRi 

^  N1/2/f  (1  +  Yx  SNRt)3/2  N  "  “ 

YlSNRl 

B2e  "  1  +  Y1  SNRt  * 


(4  .32) 


Noting  that  (1  +  02?2)  2  N  -*■  00  exp  {-5  Q  +  ■J’  SNR  )  ^  *  we  ^in<*  that 
12  v  1  i 

--C 

lim  <j>  (C)  *  e  ^  ,  i.e.,  n  is  approximately  normal  for  large  N.  With 

N+-  0 

this  important  result,  we  note  that  the  ordering  of  the  bracketed  terms 

in  the  Edgeworth  expansion  (4.20)  and  (4.21)  are  arranged  such  that 

-1/2 

successive  terms  decrease  as  increasing  integral  powers  of  N  ,  as 
shown  in  [14].  Even  with  these  important  remarks,  we  still  have  no 
feel  for  the  error  in  our  approximation.  (In  fact,  the  Edgeworth  series 
does  not  always  converge  [14].)  Consequently,  we  will  now  consider  a 
method  for  the  numerical  computation  of  H(x)  directly  from  tt'e 
characteristic  function. 

A  well-known  result  from  probability  theory  is  that  if  a  and  b 
are  continuity  points  of  a  cdf  F(.)  having  ch.f.  <}>(.)>  then 


F(b)  -  F(a)  =  lim 
T-h» 


J_  e~Jta  - 
2tt  jt 


-jta  _  -jtb 


4>(t)dt. 


(4.33) 


The  similarity  between  (4.33)  and  the  well-known  Fourier  integral  suggests 
that  a  Fast  Fourier  Transform  (FFT)  algorithm  might  be  used  to  provide 


6  H 


an  efficient  numerical  computation  of  (4.33).  With  x  =  b  and  a  -  0, 
we  may  rewrite  (4.33)  as 


F(x)  -F(0)  -lim  -r-  [ 


(4.34) 


Now,  assuming  that  p(*)  is  absolutely  integrable  and  defining 


.  .1  [  -<?(~t)  H 

T  2  it  jt 


we  have 


e’^dc, 


(4.35' 


I  =  — 

T  2tt 


i  " 

--  -  Im  [  ejxt  f(x)dx]dt 

IT  t 

0  -<*> 
oo  X 

-  i  f(x)  [  j  ■Si"-x-t-  dt]dx. 


Noting  that 


sin 

t 


-j  ,  x  <  0 


dt  **  \  0  ,  x  =  0  , 


—  x  >  0 
2  »  ~ 


we  find  that 


lim  U  =  4  -  F (0) , 
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so  that  we  may  write 


T 

.  .  1,  H  1  [  *(-t)ejxt  -  Kt)j 

F(x)  *  2  11  2 it  jt 


F(x>  *  |  ^  i 

T-k» 


)e“jxt} 


(4.36) 


Equation  (4-36)  can  also  be  found  in  [13] -  If  we  could  remove  the  Im(.) 

operation  from  inside  the  integral ,  we  could  easily  implement  an  FFT 

algorithm.  Following  NutuU  [15  ],  we  alleviate  this  problem  by  using 

"The  old  give  and  take  trick."  Choose  a(t)  such  that 

T 

,  ,  1  .<  1  f  Im  {aWe-31*)  ^ 

i)  Fa(x)  *  j  -  I1"  "  - —  dt 

T-*» 


e  L  (0,  “). 


For  such  an  a(.),  we  have 


1  r  ($(t)  - 

F(x)  =  F  (x)  -  lim  7  Im  {  **  t 

^  * 

1°  .2 


aCt))e  3-X-^dt).  (4.37) 


In  particular,  choosing  a(t)  -  e  2  ,  we  have  F&(x)  “  VM  and 


H(x)  ■  -lim 


1 1„  {  f  !  a(t)>! 

IT  C 


(4.38) 


Defining 


M  (<|>(t  )  “  a(t  ))e"jxtm 

n  m  m _ 

g(x)  “  At  l  Z 


we  may  approximate  H(x)  as 


h(x£)  =“7  Im  8(V’ 


(4.39) 


where 


,  for  t  =  (m-l)AT,  m  -  1,  2,  ...,  M,  and  x^ 


2ir  (i  -  1) 
MAt 


a  =  l,  2,  ..., 


M,  we  have 
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M 

At  l 


m=l 


(♦(t)  -  a(t  ) ) 
m  m 


t 

m 


exp{ 


-J2tt(£  -  1)  (a  -  1). 

M  ' 


(4.  AO) 


which  is  precisely  in  the  form  for  which  the  FFT  is  designed.  We  now 

need  only  to  choose  a  At  and  an  M  =  T/At  so  that  the  approximation 

(4.39)  is  a  good  one.  First,  we  note  that  g(xc)  is  periodic  with 
2  ^  2  ^ 

period  — ,  i.e.,  g(x^  +  — )  =  g(x^).  In  particular,  we  have, with 
Ax  «  2ir/T: 


'8^XM  +  m  +  l^*  m  2  ^  2  1  +1>  •••»  _1 


g (mAx) 


M, 


g(x^_j_^v,  m  —  0,  1,  2,  ...,  [ 2 1  > 


(4.41) 


t  i  n 

so  that  we  have  obtained  an  approximation  of  H(x)  for  |x|  £  — Ax. 

The  approximation  of  (4.39)  and  (4.40)  involves  errors  due  to 
truncating  the  integration  limit  as  well  as  a  sampling  error  due  to 
approximating  the  truncated  integral  with  a  Riemann  sum.  Standard 
numerical  analysis  techniques  could,  of  course,  be  applied  to 
approximate  and/or  bound  these  errors  Alternatively,  in  view  of  the 
fact  that  we  may  rf write  (4.38)  at» 


H(x) 


„  -L 


(<fr(t)  -  a (t)  )e 


-Jxt 


-dt, 


(4.42) 


we  find  that  we  are  simply  approximating  a  Fourier  integral  with  the 
FFT.  Consequently,  the  results  of  Cooley,  Lewis,  and  Welch  [16],  for 
example,  can  be  exploited  to  provide  appropriate  error  bounds. 
Following  [lb],  and  defining: 

d(t)  =  !it2_^_a£tl,  (4.43) 

00 

d  (t)  -  l  d(t  +  kT), 

P  k«-® 


and 
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00 

H  (x)  -  l  H(x  +  kX), 
P  k=-°° 


it  is  easily  shown  that 


and 


H(X£)  =  i 


d(t  )  =  "j 
ra 


dp(t)  exp{-  j2lT(£,T-  -^dt, 


0 

X 


„  /  N  j (m  -  l)uAt , 

H  (u)eJV  du. 

P 


M 


Furthermore, 

{d  (t  )}M  .  and  {H  (x  ) } .  .  are  related  by: 
p  m  m=l  P  *■ 

lAt  E  .  .  ,-12ir(m-l)U  -  1), 

Hp<V  ‘  U  l  dp(tm)eXp{  1  M  - >’ 

m“l 

and 

£  .  .  ,  i2ir(m  -  1)U  -  1), 

b  (t  *  -jAx  l  H  (x  )exp{J - * - jp* 

P  m  P  *■  r 

As  noted  previously,  we  are  making  the  approximation 

,)  =  H(xff)  =  ^  l  d<VeXp{'J27T(m  M  im 


-  1 


H(xc 


m“l 


We  will  bound  the  error  in  (4.46)  via  the  inequality: 


|H(x4)-H(xt)l  1  lH(xjt)“Hp<xt)l  +  |Hp(xJl)  -fi(xt)| 
Noting  that 

lyv^M  ^  jx  ivv*dM- 

appropriate  bounds  for  |H(x^)  -  Hp(x^)|  and  ^(t^)- 
found  and  used  in  (4.47)-  The  well-known  Markov  inequality 


(4.43) 


(4.44) 


(4.45) 


L).  (4.46) 

(4.47) 

(4.48) 

|  will  be 
states  that 


(4.49) 
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Note  that  for  v  an  even  positive  integer  the  bound  in  (4.49) 
involves  the  vth  moment  of  n»  and  that  the  bound  decreases  more 
rapidly  with  increasing  e  for  larger  v.  Since  y^Q  is  available 
from  (4.18),  the  following  bound  (4.50)  may  be  computed. 

ptNieiiV'10  (4'50) 

Using  (4.50)  and  an  inequality  given  in  [17,  p.  39],  we  find  that 


P10  .  1  -  4  e2 


|H(x)|  1  ~ fo  +  e  6  2  ’ 


(4.51) 


for  all  x  such  that  x  >  e. 


Noting  that  for  all  x  such  that  jx|  we  have  |x  +  kx|  2  X, 

and  making  use  of  (4.45)  and  (4.51)  we  find  that 

oo 

|H(x)  -H  <x)|  <  l  |H(x  -  kX)|  +  |H(x  +  kX)| 


k-1 


(for  X^2) 


4x2 


3X 


<  2p10{(jX)"10  +  (|x)”10}  +  k  e  8  e  8 

,2k-  1  VN-10 


9  2 
~X 


+  2w,„  l  <^px)-lu  +  ..j3exp{-i(^X)2}, 


‘10  .  %  '  2 
k*3 


The  series  are  easily  bounded  as 


and 


2k -1.-10  .  1  f  10,..  _  I/3V,-10 


l  (^X) 


k=3 


<  — 
-  X 


2X 

2X 


J-U,  _  ±. /jlv  'l 

x  dx  *  * 


1  2 


?  ,  1,2k -1  v,  2,  .  1  f  _2X  , 

l  exn{-^( — X)  '  —  x  I  e  °x 


k=3 


X  j 


3.x 

2X 

-h2 

s  2  8 

i— 2 e 

3X 
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We  have: 


|H(*)-H  (x)|  i2-10[(iX)"10+|<|X)'l°) 


-I  X2 


+  K§X+(i+s)e8  !• 


(4.52)  • 


,M 


which  would  be  the  only  bound  necessary  if  id  (tm) the 

M 

"pre-aliased"  version  [16]  of  (d (tm) >m=1  were  available.  Similarly, 

2k.  -  1 

noting  that  for  all  t  such  that  |t|  £  T/2,  we  have  |t  +  kl|  1  —  T, 

from  (4.43)  wo  obtain 

,  v  -  kT't  I  +  I  a  ( t  -  kT) i  +  >6(t  +  kT)  [  +  ,a(t  +  kT)  | 

~  d(t)  I  1  i 


|dp(t) 


k=l 


2k  -  1 


Now,  since  |$(t)|  and  la(c)l  are  strictly  decreasing  with  |t|, 

have 

Dl  +  laC^f^T)! 


we 


| d  ( t )  -  d(t) |  1  2  l 
p  k=l 


2k-  1 


i  t  ( i 4. (if)  |  ♦ |n<|T)  |  +  i| ♦<§T)  |  +  i| a (|t)  | ) 


0 

+  T 


f  i<?(u)  1 


du  +  | 
u  i 


I' w3”)| 

a(u)  I 


du. 


(4.53) 


2t 

2 


IT 


Recall  that  for  the  problem  at  hand,  $(t)  and  a(t)  are  given  by: 

N  r.«„ 3.3  1„  02^ 


0  C t )  =  (1  +  2  t  )  exp 


2  2,2  r33lBt  -2S26t 


2  2 

1  +  8  t 


(4.54) 


and 


a(t)  =  e 


.  w — * 
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We  have: 


dt  < 


(for  N  >  2) 


u_1  (1  +  u2)  2  exK 


"2  *2  " 


8T 


»iu 


1  +  u 


r 


1  2  2  >» 

*2  ^2P  T  f  -1, 


=  exp- 


2  2  f 

1  +  3  T  | 

J 

r 

■s 

1„ 

2  2 

--3 

,3  T 

L  \ 

2  2 
1  +  3  T 


BT 
00 

if 

2  i  2  2 
1  +  3  T 


r 

1  2  2 
,  "3„3ZT 

I  ex,-< - 


2  2r 

l1  +  l!TJ 

1  2  2 
j  n  r\  **  rr>*--\ 

exp^ 


232T2 


2  2 
i  )  du 


(x-l)_1x  2  dx 


x  dx 


2  2 
1  +  3  T 


1  2  2 

f-lS2*  T  1 

’ll t «2T2  J 


N 

7  ?  2  2  2 

(N-2)B  T  (1  +  B  T  ) 


-  I 


(4.33) 


and 


-  b^T), 


’  00  _if^ 

iiMJ  dt  <  i  [ 

t  -  T  J 


_i  t2 

dt  <  -2  e  2  =  b2(T). 

T 


(4.56) 


T  T 

Combining  (4.47),  (4.48),  and  (4.52)  through  (4.56),  we  finally  obtain: 

|H(x„)-i<xt)|  I2,10[(ix)-10+|(f«-10] 

.is2  -5  X2 

1  r  «'  .  ,1  .  1  X  ,  3  , 

+  x lo  +  +  tx)  1  1 


IlX  ,  +  |a<iT>!  +T!*<|T)|+i!a(fr)  ; 


+  WiDl(fT)  +b2(fT>1- 


(4  37) 


7.. 


Note  that  X  =  2tiM/T,  and  MAt  =  T.  For  any  range  yielding  suitably 
small  errors,  we  have: 


i<i  ^>-10  « 


J  x2 

, 1  ,  _1  ,  8 

3  +  3X)e 


<<  e 


(|x)"10 

A  x2 

8  X 


y|a(|T)i  «  h» (yT)  I 


b^l  T)  «  kC^T) 

b2(|  T)  «  |a(|T) 


Hence,  a  very  good  approximation  to  (4.47)  is 
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(4.58) 

The  approximate  bound  of  (4.58)  is  shown  in  Figure  4.4  for 
N  =  100,  CORR  =  1,  SNRi  c  (.1,  1,  10),  Ye{.l,  1)>  and  several  values  of 
M.  Since  the  bound  is  decreasing  with  increasing  N  and  decreasing 
|  CORR  |  ,  Figure  4.4  suggests  that  for  T  =  40tt,  N  100,  |C0RR|  <_  1, 

.1  <  SNR.  <  10,  .1  <  y  <  1,  and  M  =  1024,  an  error  of  less  than  10 
will  result.  For  this  choice  of  T,  we  have  Ax  =  .05,  which  should 
provide  both  acceptable  resolution  and  error  in  approximating  H(x) 
by  H(x  )  in  (4.4b).  Figure  4.5  illustrates  the  numerical  results  for 
parameter  choices  corresponding  to  those  of  Figure  4.3.  A  careful  com¬ 
parison  of  Figures  4.3  and  4.5  reveals  several  discrepancies,  indicating 
that  the  numerical  inversion  technique  with  its  associated  error  bound 


is  indeed  worthwhile. 
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Although  both  the  Edgeworth  approximation  (4.22)  and  the  numerical 
inversion  computation  (4.46)  are  useful,  their  usefulness  decreases  when 
considering  the  tails  of  the  cdf  F(x).  Consequently,  a  "tight"  bound 
on  the  tails  of  F(x)  would  be  desirable.  The  Markov  inequality  (4.59) 
can,  of  course,  be  used,  but  the  so-called  Chernoff  bound  [17,  pp.  118-122] 
decreases  exponentially,  whereas  (4.49)  decreases  only  as  e  V.  Also, 
the  distribution  under  consideration  here  is  skewed  (since  f4  0),  a 
property  which  is  lost  when  using  bounds  such  as  (4.49) 

Consider  for  X  >_  0: 

00 
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for  all  X  >  0,  such  that  «  (-j A)  exists.  From  the  development  of 


(4.8)  it  is  easily  seen 
Similarly,  for  X  0: 


that  $  (— j  A)  exists  only  for  A  8  <  1. 
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for  all  A  0  such  that  <>  (j A)  exists. 
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Defining  h(A)  =  Lot>  ^(-jA)  (A  ft  •'1),  we  have 

P{n  >  e}  <_  e  At"  +  for  0  £  A  <  1/3, 

and 

P{n  <  e}  <  e"A£  +  h(X),  for  -1/3  <  A  <_  0. 

We  now  seetc  to  minimize  the  quantity  b(A)  =  exp{h(A)  -  Ae}  with  respect 
to  A.  We  have: 


b’(A)  =  (h'(A)  -  e)  b ( A )  =  0 

iff  h' (A)  -  e,  where  h'(A)  -  e  will  result  in  the  minimum  b(A> 

if 

i)  there  exists  a  A  such  that  h' (A)  =  e 

and 

ii)  b" (A)  |  >  0  . 

e*h' (A) 
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The  bounds  of  Interest  become 
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(4.60) 


P(n  <  h' (A; }  <  eh(A)  "  Ah'(A),  for  ~  <  A  <  0. 

We  note  that  (h'(A):  0  <  A  <  |)  =  (x:  0  <  x  <  ~)  and 
{h'(A);  _i  <  \  <  o)  =  (x:  <  x  <_  0).  The  bounds  of  (4.60) 

illustrated  in  Figure  5.6  for  |h'(A)|  <_  10,  CORR  =  1, 

Ne{2,  10,  100,  1000),  SNRjd.i,  1,  10),  and  y  e  (.1,  1).  Note 
the  bounds  of  (4.60)  are  decreasing  with  |C0RR|  and  that 

P(n  <  e) |  =  P(n  >  -t) |  » 

-CORR  CORR 
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We  now  have  a  relatively  complete  characterization  of  the  cdf  F^GO* 


z  -  V. 


where  n  is  the  standardized  r.v.  n 


In  order  to  use 


these  results  for  the  r.v.  z,  we  note  that 


/*  '  U: 
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(4.62) 
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V  Conclusion 

In  this  report,  we  have  studied  two  predistorted-replica 
correlation  receivers:  The  FLAC  receiver  and  the  ARSAC  receiver.  A 
Doppler  analysis  was  performed  for  each,  the  results  of  which  indicate 
that  for  signal  time  bandwidth  products  TW  such  that 


the  difference  in  sensitivity  is  not  dramatic;  however,  for  larger 
time-bandwidth  products  the  ARSAC  receiver  vastly  outperforms  the 
FLAC  receiver. 

Sections  II-C  and  III-C  pointed  out  that  the  FLAC  concept  can 
perform  satisfactorily  in  a  multiple  user  environment;  whereas  the 
ARSAC  concept  will  not. 

Sections  II-B  and  III-B  dealt  with  discrete  time  structures  for 
the  FLAC  and  the  ARSAC  systems,  respectively.  In  fact,  a  digital 
realization  of  the  FLAC  receiver  can  be  made  requiring  only  two 
additions  and  one  multiply  per  output  sample.  The  output,  y  ,  for  the 
FLAC  receiver  is,  from  Section  II-B: 


or 


L-l 
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(5.1) 


Rn-L+l(Rn  +  l  "  Rn-2L  +  1*  +  V 


(5.2) 


indicating  that  if  y^  “  0.  and  Rn  =  0  for  all  n  <  £,  (5.2)  is 
equivalent  to  (5.1)  for  all  n  =  £,  H+l,  £.+2,.... 
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It  is  to  be  emphasized  that  the  assumptions  made  about  the  channel 
imply  that  there  is  no  intersymbol  interference,  an  assumption  which 
is  certainly  restrictive.  For  the  known  signal  case,  seve"al  authors 
have  shown  that  the  optimal  structure  is  a  matched  filter  followed  by 
a  transversal  filter  [3].  It  is  conjectured  here  that  similar  per¬ 
formance  improvement  can  be  obtained  by  following  the  predistorted 
replica  correlation  receiver  by  a  transversal  filter. 

Finally,  the  noise  performance  for  the  family  of  predistorted 
replica  correlation  receivers  was  dealt  with  in  Section  IV.  The 
analysis  there  assumed  that  the  noise  component  of  the  received  process 
was  additive  Gaussian  noise,  such  that  the  samples  of  the  received 
process  were  uncorrclated. 
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